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Context. Distances to most star forming regions are determined using kinematics, through the assumption that the observed radial 
velocity arises from the motion of the source with respect to the Sun resulting from the differential rotation of Galaxy. The primary 
challenge associated with the application of this technique in the inner Galaxy is the kinematic distance ambiguity. 
Aims. In this work, we aim to resolve the kinematic distance ambiguity towards a sample of 6.7 GHz methanol masers, which are 
signposts of the early stages of massive star formation. 

Methods. We measured 21 cm Hi absorption spectra using the Very Large Array in C and CnB configurations. A comparison of the 
maximum velocity of H i absorption with the source velocity and tangent point velocity was used to resolve the kinematic distance 
ambiguity. 

Results. We resolved the distance ambiguity towards 41 sources. Distance determinations that are in conflict with previous measure- 
ments are discussed. The NE2001 spiral arm model is broadly consistent with the locations of the star forming complexes. We find 
that the use of vertical scale height arguments to resolve the distance ambiguity can lead to erroneous classifications for a significant 
fraction of sources. 

Key words. Galaxy: kinematics and dynamics - Galaxy: structure - Hn regions - Radio continuum: ISM - Radio lines: ISM 



1. Introduction 

Measuring distances towards massive star forming regions is im- 
portant both for characterizing the physical properties of these 
regions, and for studying the spiral structure of the Galaxy. 
Trigonometric parallax is the most direct and model indepen- 
dent technique for measuring distances to astronomical objects. 
. £^ | However, it is difficult to measure trigonometric parallax towards 
. a large sample of massive star forming regions due to their large 
' distances. The popular method for deriving the distance to a 
source arbitrarily located in the disk of the Milky Way is to as- 
sume that the measured radial velocity of the source arises from 
its differential rotation in the Galax y. Then, by usi ng a model 
for the rotation of the Galaxy (e.g. IClemenslll98"5t C85 here- 
after), one can obtain an estimate of the kinematic distance to 
the source. One of the main problems associated with kinematic 
distances is that the distance solution is double valued in the 
inner Galaxy. Thus, a given radial velocity for a source in the in- 
ner Galaxy corresponds to two distances - a near distance and a 
far distance. The two solutions are equidistant about the tangent 
point, which is the location where the line of sight is a tangent to 
the assumed circular orbit of a source about the Galactic center. 
This problem is called the kinematic distance ambiguity (KDA) 
and is the primary obstacle towards a straightforward application 
of kinematic distances to study Galactic structure. 

One of the techniques to resolve the KDA is to measure an 
absorption spectrum towards the source of interest. Since only 
material in front of the target can give rise to absorption, one 
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can use the velocities of the absorption features t o distinguish 
betwe en the near and far distance (see Fig. 1 of iKolpak et al.l 
2003 for a schematic). Thus, if absorption features are seen at 
velocities between that of the source and of the tangent point, 
the source has to be at the far distance since the clouds produc- 
ing absorption at these velocities will be located between the 
near and far distance solution points (and will consequently be 
further from the Earth than the near distance point). For the tech- 
nique to work, the tracer used to measure the absorption has to 
be abundant in the Gala xy. The two popular tracers used for this 
purpose are H i at 2 1 cm ( Kuc har & Baniall994tlFish et al.l2003k 
I Kolpak et alj|2003l). and formaldehyde (CH?Q) at 6 cm dWilsonl 
ll972tlArava et al .1120021: TWatson et all2003USewilo et alj2004l) . 

Each of the above tracers has its advantages and disadvan- 
tages. Formaldehyde is less widespread than Hi, and this can 
lead to an incorrect resolution of the KDA. One way this can 
come about is simply that there is a low formaldehyde abundance 
in clouds between the near and far distance points. The resulting 
non-detection evidently leads to incorrect classification of a far 
distance source to be at the near distance. The low excitation 
temperature of the 6 cm formaldehyde transition means that no 
discrete background source is necessary, as formaldehyde ab- 
sorption lines can be detected against the Cosmic Microwave 
Background (CMB). In contrast, H i is much more abundant and 
widely distributed in the interstellar medium, so that a cloud 
lacking HI is a rare occurrence. Consequently, any KDA reso- 
lution that is obtained in this way is likely to be correct. It is 
true that absorption spectra from single dish telescope observa- 
tions are difficult to interpret due to emission from gas spatially 
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extended be yond the boundaries of the background continuum 
source (e.g. iDickev & LockmanllT990i) . This confusion can be 
mitigated through interferometer observations, as discussed in 
Section 3. Another potential is that young Hn regions may lack 
21 cm continuum emission against which to measure H i absorp- 
tion. In some cases, where the emission is optically thick, the 
continuum may be more easily detectable at 6 cm compared to 
21 cm, but at least for observations of H i, lack of sufficient con- 
tinuum will be evident. Thus, the consequence is simply a lack 
of data from such sources. 

Masing transitions of molecules such as methanol, OH and 
water are frequently associated with early phases of massive star 
formation. Methanol masers at 6.7 GHz, the brightest of Class 
II methanol masers, are unique in that they appear to be exclu- 
sively associated with massive star formation (e.g. Mini er et"ai] 
2003; iPandian & Go ldsmith 2007). Thus, methanol masers are 
potentially useful to further our understanding of Galactic spi- 
ral structure. This is bolstered by the recent measurement of the 
distance to 12.2 GHz methanol masers (which are the second 
brightest C lass II methan ol masers) in W30H, and the Perseus 
spiral arm dXu et alj|2006l) . It is however, not possible to apply 
the absorption spectrum technique to resolve the KDA towards 
a majority of methanol masers in the inner Galaxy. This is be- 
cause the sources immediately associated with the masers of- 
ten do not have detectable radio continuum (especially at long 
wavelengths) due to the young age of the massive young stel- 
lar objects (MYSOs). For example, only about 10% of the 6.7 
GHz methanol masers detected in the blind Arecibo survey hav e 
detectable continuum at 21 cm ( Pandian & Goldsmith 2007). 
However, it is well known that massive stars form in clustered 
environments. Thus, while the number of methanol masers di- 
rectly associated with radio continuum may be small, it is likely 
that a larger fraction will be in a cluster hosting a more evolved 
MYSO that has associated radio continuum emission. Here, we 
report on the results of Hi observations with the Very Large 
Array (VLA) to resolve the KDA towards a sample of 6.7 GHz 
methanol masers that either have direct association with 21 cm 
continuum, or are located in a cluster hosting a 2 1 cm continuum 
source. 

2. Sample Selection 

The methanol maser sampl e was selected from the General 
Catalog of methanol masers (Pest alozzi et al.ll2005l) . We first re- 
stricted ourselves to methanol masers in the first Galactic quad- 
rant with Galactic longitude I > 10°. At low Galactic longitudes, 
a large fraction of sources in the line of sight have a signif- 
icant component of their circular motion (around the Galactic 
center) oriented tangential to the Sun. Consequently, for these 
sources, radial peculiar motions give rise to large errors in the 
kinematic distance. Sources near the tangent point do not have 
this problem, but the difference between the near and far distance 
is not as large compared to sources at larger Galactic longitudes. 
Furthermore, these sources are at low declination where the sen- 
sitivity of the VLA is poorer. 

We searched for point sources i n the NRAO VLA Sky 
Survey (NVSS; ICondon et ail Il998l) within ~ 2' of each 
methanol maser. We placed a cut-off of 30 mJy for the flux of 
the continuum source in the interest of observing time required 
to obtain an absorption spectrum with good signal to noise ratio. 
We then looked for radio recombination line (RR L) associations 
with the continuum sources (using the catalog o flLock man 1989; 
L89 hereafter) to verify their nature as Hn regions. The veloc- 
ities of the recombination lines were also used to check the as- 



sociation of the H ii region with the cluster hosting the methanol 
maser. Sources with no RRL information in L89 were observed 
with the National Radio Astronomy Observatory'^] Green Bank 
Telescope (GBT) for this purpose. Sources with no RRLs were 
discarded as they are likely to be continuum sources other than 
Hn regions that are projected against our target masers. These 
selection criteria gave a sample of 64 methanol masers that we 
observed with the VLA. 



3. Observations 

3.1. GBT Observations and Data Reduction 

The observations were made in June 2005 using the C-Band re- 
ceiver of GBT. The GBT spectrometer was used in 16 interme- 
diate frequency (IF) mode to cover 8 frequencies in two polar- 
izations. The IFs were centered on the hydrogen recombination 
lines H103q' to HllOff. The frequencies adopted for the tran- 
sitions were 5931.545 MHz, 5762.881 MHz, 5600.551 MHz, 
5444.261 MHz, 5293.733 MHz, 5148.703 MHz, 5008.923 MHz 
and 4874.158 MHz respectively for the HlOSc to HI 10a lines. 
Each band had a bandwidth of 12.5 MHz with 4096 channels, 
giving a total velocity coverage of ~ 650 km s _1 with ~ 0.2 
km s _1 resolution. The observations were made in frequency 
switched mode with a 3 MHz frequency throw. The integration 
time per source was typically 15 minutes, with shorter integra- 
tions (5-10 minutes) being used for sources with strong contin- 
uum. The system temperature was usually around 25 K. A point- 
ing and focus check was done at the beginning of each observa- 
tion. The flux scale was set using noise diodes in the receiver. 
Since the intent of the observations was to detect RRLs and their 
velocities, we did not do any additional flux calibration using a 
calibrator. We estimate the flux scale to be accurate to around 
10%. 

The data were reduced using GBTIDL. Frequency switched 
observations yield two sets of difference spectra for each IF. 
Each difference spectrum was baseline subtracted using a 3rd 
order polynomial. The baseline subtracted difference spectra of 
the IFs that were not affected by radio frequency interference 
were then summed together, taking into account the different ve- 
locity resolutions at the different frequencies. The resulting dif- 
ference spectra were then folded to give the final spectrum of the 
source. This technique gives spectra with fairly good baselines 
that have a signal to noise ratio as much as a factor of 4 higher 
than what would be achieved using a standard position switched 
observation of a single recombination line. A Gaussian fit was 
performed on the detected RRLs to determine their velocity with 
respect to the local standard of rest (LSR). 

3.2. VLA Observations and Data Reduction 

A total of 63 different fields (covering 64 sources) was observed 
with the VLA between June and October 2005. Sources at de- 
clinations below -15° were observed in the CnB configuration, 
while the other sources were observed in the C configuration. 
The quasar 3C 286 was used as a primary flux calibrator. The 
correlator was configured to have a bandwidth of 1 .5625 MHz 
with 256 channels in single IF mode. The calibrators were ob- 
served at +250 km s _1 and -180 km s _1 (with respect to the 
rest frequency of H i in the LSR frame of reference, assumed to 
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be 1420.40575 MHz) to avoid contamination from Galactic Hi, 
while the targets were observed at a velocity of +50 km s _1 . The 
total integration time per field was 10-25 minutes depending on 
the NVSS flux of the continuum source. 

The data were reduced using the Astronomical Image 
Processing System (AIPS) of NRAO employing standard pro- 
cedures. The initial calibration was done using the "channel 0" 
data, the solutions from which were copied to the line data. For 
each source, we summed line free channels to create continuum 
datasets. We then CLEANed and self-calibrated the continuum 
data, and copied the resulting solutions to the line data. We were 
conce rned about CLEANing the line data since iKolpak et al.l 
(2003) reported that the procedure introduced spurious absorp- 
tion features from the emission features that are not resolved by 
the interferometer, leading them to extract spectra from the dirty 
image cubes. While the latter procedure works for strong con- 
tinuum sources that are not affected by the sidelobes of other 
sources in the field, it fails for the weaker sources in our sample. 

Hence, we did an extensive test of CLEANing the line data. 
We used the strong sources as a reference to check for any intro- 
duction of spurious absorption features. We found that as long 
as the CLEAN is restricted to ~ 2<x of the noise in a channel 
image, the procedure works without any artifacts being intro- 
duced. We also found that using boxes to restrict the area that 
is CLEANed did not work well, as it introduced Hi emission 
from outside the CLEAN boxes (which are centered around the 
continuum source) into the targets. Although the emission fea- 
tures lie outside the CLEAN boxes, they enter the CLEAN boxes 
through the sidelobes of the dirty beam. The use of CLEAN 
boxes would thus force the task to clean the sidelobe emission 
instead of the peak which lies outside the CLEAN boxes. The 
exception is when one is dealing with fields that have both weak 
and strong point sources that are embedded in a large extended 
structure. In these cases, a full field CLEAN did not reproduce 
the weak point sources, and introducing boxes to restrict the area 
CLEANed gave much better results. We also found that remov- 
ing the continuum prior to the CLEAN (using the task UVLSF) 
yielded better results compared to CLEANing the fields with no 
continuum subtraction. All the spectra presented in this paper are 
derived from CLEANed image cubes. The CLEANed data were 
then divided by the continuum image to derive optical depth in- 
formation. 

4. KDA Resolution and Distance Calculation 

To resolve t he KD A, we used the approach followed by 
Kolp ak et al.l d2003). who performed a simulation to explore the 
relation between the quality of KDA resolution and the tan- 
gent point velocity, v,, source velocity, v s , and the highest ve- 
locity of a H i absorption feature, v a , taking into account ran- 
dom velocities, streaming motions and cloud-cloud dispersions. 
They found that the KDA can be resolved with good confidence 
unless (v, - v„) < (v, - v a ) + 15 km s and v t — v a < 20 
km s -1 (these are sources with velocities close to the tangent 
point velocity). Further, within the ambiguous region, the line 
(v; - v a ) = (v ( - v,)/2 separated the near distance an d far distance 
sourc es at the confidence level of 50% (Fig. 2 of IKolpak et al.l 
120031) . We use the same criteria, and assign sources to the near 
or far distance with the quality grade of A and B depending on 
whether the source lies outside or within the ambiguous region, 
respectively. 

To calculate distances, we use two approaches. The first is 
the classical approach, employing the source velocity, v s , to cal- 
culate the distance using the C85 rotation curve with (Rq, ©o) = 



(8.5, 220). We prefer to use the C85 rotation curve as all our 
sources are in the first quadrant, and the C85 curve takes into ac- 
count peculiar motions in this region, s uch as the dip in velocit y 
near a Galactocentric radius of 3 kpc (Bur ton & GordorJll978l) . 
The errors in the distance are calculated using the distances cor- 
responding to v s + 10 km s to allow for a the velocity of the 
region differing from that arising from pure circular rotation. 

However, parallax measurements of massive star forming re- 
gions show a consistent bias in that kinem atic distanc es over- 
predict the true distance to the source rfReldl l2007h . Using 
VLBI measurements of the proper motion of Sagittarius A*, the 
value of @o/^o has been measured to be 29.5 km s _I kpc -1 
with 0.3% accuracy, assu ming that Sagittarius A* is stationary 
dReid & Brunthalerll2004 . This translates to a O of 250.8 km 
s _1 for an Rq of 8.5 kpc. The use of these new parameters still 
yields significant systematic peculiar motions for the young mas- 
sive star forming regions as seen in the frame of reference that is 
rotating with the G alaxy . Doing a global fit to reduce the pecu- 
liar motions. iReidl ((2007) find that massive star forming regions 
are born rotating ~ 13 km s slower than the Galactic rotation 
curve, and that the motion o f the Sun towards the LS R is con- 
sistent with that derived by iDehnen & Binnevl d 19981) from the 
Hipparcos data. 

Hence, in an attempt to make kinematic distances closer 
to the parallax distance, we used the following approach: we 
first calculated the heliocentric velocity from the current LSR 
model (which assumes a solar motion of 20 km s towards 
(a, 6) = (19'\ 30°) in 19 00 coordinates). We then recalculated 
an LSR velocity using the lDehnen & Binnevl (I 1998b model with 
(U Q , Vq, Wo) = (10.0, 5.25, 7.17). Using the new velocity (typ- 
ically less than the old velocity by ~ 10 km s~'in the first and 
second Galactic quadrants), we calculated the distance using the 
C85 curve with (R , O ) = (8.5 kpc, 250.8 - 13 km s 1 ). To 
calculate uncertainties, instead of using the traditional approach 
of calculating kinematic distances corresponding to v s + 10, we 
used a ±10 km s variation to the component of the source ve- 
locity that is radial to the Galactic center. This is because the 
peculiar motions of massive star forming regions, in addition to 
being ~ 13 km s slower to Galactic rotation, are preferentially 
oriented in the dire ction of the Ga lactic Center, as is seen for 
example in W30H dXu et alj|2006l) . We also included a mini- 
mum uncertainty of ±5 km s _1 to the source velocity to take 
into account random motions, since the above prescription will 
not yield any uncertainty to the LSR source velocity when the 
source is near the tangent point. 

The difference between the distances calculated using the 
two methods is typically small (~ 0.2 kpc), but can be larger 
than 0.5 kpc, especially at velocities near the tangent point ve- 
locity, and in the second Galactic quadrant. 

5. Results 

The continuum sources associated with (either directly or 
through a cluster) 36 methanol masers were detected with ad- 
equate signal to noise ratio to resolve the KDA. The remaining 
28 sources were typically too resolved to extract spectra with 
good enough signal to noise ratio. Treating multiple recombina- 
tion line velocities (towards the same line of sight) as separate 
sources, we have a total of 41 sources. The properties of these 
sources are listed in Table 1. Column (1) contains the source 
names. Sources which have two recombination line velocities 
in the same line of sight are suffixed A and B for distinction. 
Column (2) lists the NVSS source that is associated with the 
methanol maser. Column (3) shows the velocity extent of the 
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maser emission, and column (4) shows the systemic velocity, 
v s derived from the radio recombination lines. Column (5) lists 
the KDA resolution specifying whether the source is at the near 
distance (N), far distance (F), tangent point (T), or whether the 
KDA is unresolved (?). Column (6) shows the distance d, calcu- 
lated from the conventional method using the C85 rotation curve, 
while column (7) shows d' , which is the modified kinematic dis- 
tance calculated using the methodology described in section 4. 
Column (8) indicates the quality of the KDA resolution (A or B), 
and column (9) lists the rms noise in the 21 cm continuum im- 
age. Errors for d are not listed when v s + 10 exceeds the tangent 
point velocity, v f , since such velocities are forbidden, and cannot 
be converted to a meaningful distance. 

Figure 1 shows the continuum images, H i optical depth, and 
GBT recombination line data (where applicable) for our sample. 
The optical depth spectra also indicate the source velocity, v s , 
and the tangent point velocity, v t , as short and long dashed lines, 
respectively. 

5. 1 . Notes on Selected Sources 

10.32-0.15 and 10.62-0.38 - The methanol maser emission in 
10.32-0.15 and 10.62-0.38 range from 4 to 17 km s 1 and -10 
to 7 km s _1 respectively. The recombination line velocities from 
L89 for the two sources are 11.0 and 0.4 km s _1 respectively. 
Both sources are associated with the W3 1 complex, the distance 
to which has been a matter of contention. Absorption line exper- 
iments have not discovered any absorption features beyond ~ 45 
km s . A purely straightforward resolution of the KDA would 
put both sources at the far distance. However, it is intriguing that 
no H i absorption is seen over a velocity interval corresponding 
to ~ 7 kpc, leading to suggestions that the source is at the near 
distance with the systemic velocity diffe ri ng from circular mo- 
tion by ~ 35 km s" 1 (e.g. IWilsonll 19741) . ICorbel & Eikenberrvl 
(2004) suggest that the lack of absorption features in both atomic 
and molecular gas is due to a large gas hole between Galactic 
longitudes of 5°and 25° at velocities greater than 45 km s _1 . 
However, there are a number of sources in the same Galactic 
longitude range that do show H i absorption features between 45 
km s^ 1 and the tangent point velocity (as can be seen from Fig. 
1). This casts doubt on this explanation, although it is possible 
that there is a lack of H i in the specific direction of W3 1 . 

The situation is complicated by the observation of the nearby 
luminous blue variable star LBV 1 806-020 to be a t the far dis- 
tance of 15.1+1* kpc (ICorbel & Eikenberrvl 120041: as inferred 
by the presence of ammonia absorption at 73 km s _1 ) and the 
derivation of a spectrophoto metric distance to the H n region 
G10.2-0.3 to be 3.4 + 0.3 kpc dBlum et alJIiOOlh . although both 
sources have similar LSR radial velocities. This suggests a su- 
perposition of complexes at vastly different distances which co- 
incidentally have very similar radial velocities. A parallax mea- 
surement of the distances to the regions inside this complex is 
thus required to determine the true distances and the peculiar 
motions of the various sub-regions. In light of the fact that Hi 
absorption is seen at velocities higher than 45 km s at sim- 
ilar Galactic longitudes (10.47+0.02 for example), we suggest 
that both these regions are at the same distance as G10.2-0.3, 
3.4 kpc. The difference between the LSR velocity expected for a 
source at that distance, and the observed LSR velocity can then 
be attributed to peculiar motions. For example, W30H has a 
large peculiar motion of 17 km s towards the Galactic center in 
additi on to rotating 14 km s _1 slower than the Galactic rotation 
curve (IXu et aD HOOg). We stress that an independent study that 



is not based on kinematics is required to resolve the conflicting 
ideas on the true distance to the regions in the W3 1 complex. 

12.20-0.10- There are two recombination line velocities de- 
tected here at 27.0 and 53.6 km s _1 . The methanol maser emis- 
sion extends from 2 to 15 km s _1 and is thus associated with 
the 27.0 km s _1 Hn region. The continuum image shows mul- 
tiple components, which could give rise to the two recombina- 
tion line systems. An alternate possibility is that one of the Hn 
regions giving rise to the recombination lines is optically thick 
and is undetected at 21 cm. This is possible as the measurements 
of L89 were made at 3.6 cm, and an optically thick Hn region 
can be 30 times weaker at 21 cm compared to 3.6 cm. The Hi 
spectra from all components of the 21 cm image show absorp- 
tion features between the systemic velocity and the tangent point 
placing the source(s) at the far distance. Hence, while it is pos- 
sible that one of the H n regions is not detected in the continuum 
image at 21 cm, we tentatively place both regions including the 
methanol maser at the far distance. 

12.68-0.18 - The GBT spectrum (at 6 cm) shows two re- 
combination line systems, at 32.0 and 57.4 km s , the 32.0 km 
s component being weaker. The methanol maser emission ex- 
tends from 50 to 62 km s _1 and is thus associated with the 57.4 
km s^ 1 recombination line. The 21 cm continuum image shows 
a very extended structure to the source, the component marked 
'X' being closest to the methanol maser position. The H i spec- 
trum for this source shows absorption features up to ~ 55 km 
s , which is consistent with it being associated with the 57.4 
km s region and it being at the near distance. The signal to 
noise ratio at the other continuum regions is too low to deter- 
mine whether there is another type of H i spectrum that is asso- 
ciated with the 32.0 km s -1 component. Hence, we classify the 
methanol maser to be at the near distance, and do not make any 
claim about the second H n region in the line of sight. 

18.06+0.08 - The GBT spectrum shows two velocity com- 
ponents, at 16.4 and 60.1 km s . The methanol maser, whose 
emission ranges from 45 to 57 km s _1 , is associated with the 
60.1 km s component. The 21 cm continuum is dominated by 
a single source, although it is possible that some of the structure 
(especially the extension to the north-west) arises from a super- 
position of a second source. In the absence of a 6 cm image, we 
are unable to determine whether both H n regions are seen at 21 
cm. Only the central source has adequate signal to noise ratio 
to resolve the KDA, the spectrum from which is shown in Fig. 
1 . The spectrum would place the source at the far distance for 
either RRL velocity, and thus we indicate the two possibilities in 
Table 1 without making any claim regarding which component 
is seen in our image. 

23.43-0. 18- The GBT spectrum shows a prominent compo- 
nent at 104.5 km s _1 which is associated with the maser emission 
which spans a velocity range between 94 to 113 km s _1 , along 
with a much weaker component at 60.6 km s _1 . The tangent 
point velocity is 1 13.6 km s" 1 . This sour ce is classifie d at th e 
near distance by ISewilo et al.l d2004l) and iKolpak et al.1 d2003l) . 
However, our absorption spectrum shows an absorption compo- 
nent at 1 16.5 km s _1 . We thus assign the far distance with quality 
grade B to this source. 

24.50-0.03 - This source, with a systemic velocity of 109.2 
km s _1 (maser emission ranges from 108 to 116 k m s _I ) is 
classified at the near distance by Sew ilo et~ai] (12004). The tan- 
gent point velocity at this Galactic longitude is 112.9 km s _1 . 
Our H i spectrum shows an absorption feature at around 130 km 
s~'detected at greater than 3cr level. We thus classify this source 
to be at the far distance with quality grade A. 
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Table 1. 6.7 GHz Methanol Maser Properties 



Source 


NVSS 


Maser (v„„„, v„ M ,) 


v., 


KDA 


d 


d' 


Quality 


a 


Name 


Source 


(km s" 1 ) 


(kms- 1 ) 


Resolution 


(kpc) 


(kpc) 


Assignment 


mjy beam -1 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


10.32-0.15^ 


180900-200456 


(4, 17) 


11.0 


N 


3.4 + 0.3 


3.4 + 0.3 




8.0 


10.47+0.02" 


180835-195219 


(58, 77) 


70.1 


F 


10.6 ± +0.4 


10.8 + 0.3 


A 


5.3 


10.62-0.38" 


181028-195548 


(-10, 7) 


0.4 


N 


3.4 ±0.3 


3.4 ± 0.3 




1.5 


10.95+0.02 


180942-192629 


(23, 25) 


18.5 


F 


14 2 + ' 3 


14.3 ± 1.2 


A 


1.3 


11.03+0.06 


180939-192116 


(15,21) 


12.6 


F 


15.0+1.3 


15.0+ 1.3 


A 


1.3 


11.93-0.61 


181401-185324 


(30, 44) 


39.3 


N 


4.2 + 0.8 


4.0 ± 0.7 


A 


2.2 


11.99-0.27 


181245-184039 


(60, 61) 


58.2 


F 


11.3 + 0.6 


11.5 + 0.6 


A 


1.4 


12.20-0. 10A 


181242-182417 


(2, 15) 


27.0 


F 


13.5 ± 1.0 


13 6 +l 

"■"-0 7 


A 


1.3 


12.20-0. 10B 


181242-182417 




53.6 


F 


11.6 + 0.6 


11.8 + 0.6 


A 


1.3 


12.68-0.18 


181355-180124 


(50, 62) 


57.4 


N 


5.1+0.6 


4.9 ± 0.6 


A 


5.5 


12.79-0.19 


181414-175543 


(31,41) 


35.8 


N 


1 0+0.6 


3.6 ± 0.7 


A 


2.5 


12.90-0.26 


181444-175223 


(35, 47) 


28.6 


F 


13.4+1.0 


13.6 + 0.9 


A 


2.5 


13.18+0.06 


181405-172832 


(48, 49) 


54.1 


N 


4.8 + 0.6 


4.6 + 0.5 


A 


2.5 


14.09+0.10 


181545-163858 


(4, 17) 


16.3 


F 


14.7 ± 1.1 


14.8 ± 0.9 


A 


1.5 


18.06+0.08A^ 


182335-130920 




16.4 


F? 


14.7 + 0.9 


14.8 ± 0.8 


A 


2.0 


18.06+0.08B*" 


182335-130920 


(45, 57) 


60.1 


F? 


11.7 + 0.5 


11.9 + 0.4 


A 


2.0 


18.46-0.00 


182436-125106 


(47, 50) 


56.5 


F 


11.9 + 0.5 


12.1+0.4 


A 


2.2 


19.48+0.15 


182604-115234 


(20, 27) 


19.8 


F 


14.4 ± 0.8 


14.5 + 0.7 


A 


2.2 


19.61-0.13 


182716-115322 


(52, 54) 


58.6 


F 


11.8 + 0.5 


12.0 + 0.4 


A 


1.7 


23.43-0.18 


183438-083043 


(94, 113) 


104.5 


F 


9.2 


9.6 + 0.3 


B 


3.3 


24.50-0.03 


183606-073106 


(108, 116) 


109.2 


F 


8.6 


q -2+03 
y - J -0.5 


A 


1.5 


24.67-0.14 


183652-072453 


(111, 114) 


111.4 


N 


7.2 


6.4 ± 0.5 


B 


2.4 


24.78+0.08 


183610-071124 


(106, 115) 


108.6 


N 


6.8 


6.1+0.4 


B 


1.6 


25.70+0.04A 


183804-062523 




53.3 


F 


11.8 + 0.6 


12.0 + 0.4 


A 


1.0 


25.70+0.04B 


? 


(89, 101) 


102.0 


? 








1.0 


27.28+0.15 


183610-071124 


(34, 37) 


36.3 


F 


12.7 + 0.7 


12.9 + 0.5 


A 


0.8 


28.30-0.38Aj. 


184415-041757 




44.8 


F 


12.1+0.6 


12.3 + 0.5 


A 


0.8 


28.30-0.38B~ 


184422-041746 


(79, 93) 


75.6 


F 


10.3 ± 0.6 


10.6 ± 0.3 


A 


0.8 


30.70-0.07 


184735-020143 


(85, 90) 


101.0 


F 


8.5 


9.0 ± 0.4 


A 


4.4 


30.76-0.05 


184738-015731 


(88, 95) 


91.6 


F 


8.5 


9.0 ± 0.4 


B 


4.4 


30.82-0.05 


184746-015451 


(91, 110) 


105.7 


F 


8.5 


9.0 ± 0.4 


B 


4.4 


32.97+0.04 


185124+000412 


(89, 93) 


93.2 


F 


8.3 


8.7 ± 0.4 


B 


2.0 


33.13-0.09 


185207+000819 


(71, 81) 


93.8 


F 


8.2 


8.7 ± 0.4 


B 


1.3 


33.40+0.01A 


185220+002552 




72.3 


F 


9.8 ± 0.7 


10.0 + 0.4 


A 


1.3 


33.40+0.01B 


185216+002531? 


(96, 107) 


105.4 


T 


7.1 


7.1 


A 


1.3 


41.10-0.22 


190711+070956 


(63, 64) 


59.4 


F 


9.1+0.8 


9.4 ± 0.4 


A 


1.9 


43.87-0.77 


191426+092236 


(47, 53) 


55.0 


F 


8.8 + 0.8 


9.1+0.5 


A 


2.7 


45.07+0.13 


191320+105054 


(57, 60) 


63.0 


F 


7.7 


8.0 + 0.5 


B 


1.4 


48.90-0.27^ 


192214+140319 


(63, 73) 


66.5 


? 


7.0+1.5 


7.0+ 1.5 




1.4 


50.29+0.69 


192127+154406 


(26, 33) 


25.9 


F 


9.3 + 0.7 


9.3 + 0.3 


A 


1.3 


81.87+0.78" 


203837+423759 


(0, 13) 


1.7 


F 


2.0 


1.2 


B 


2.9 



Column (1) contains the source names. Sources which have two recombination line velocities in the same line of sight are suffixed A and B for 
distinction. Column (2) lists the NVSS source that is associated with the methanol maser. Column (3) shows the velocity extent of the maser 
emission, and column (4) shows the systemic velocity, v s derived from the radio recombination line. Column (5) lists the KDA resolution 
specifying whether the source is at the near distance (N), far distance (F), tangent point (T), or whether the KDA is unresolved (?). Column (6) 
shows the distance d, calculated from the conventional method using the C85 rotation curve, while column (7) shows d', which is the modified 
kinematic distance calculated using the methodology described in section 4. Column (8) indicates the quality of the KDA resolution (A or B), and 
column (9) lists the rms noise in the 21 cm continuum image. Errors for d are not listed when v, + 10 exceeds the tangent point velocity, v,, since 
such velocities are forbidden, and cannot be converted to a meaningful distance. 

a Distance from lBlum et all d200ll) 

* Source has ambiguity between two recombination lines 

The continuum source associated with the A and B sources may be reversed 

Distance from lGenzel et alj dl981l) 
' d' for this source is the distance to the tangent point 



24.67-0.14 - The systemic velocity of this source is 111.4 
km s , while the tangent point is at 112.9 km s . The highest 
velocity of an absorption feature is 110.5 km s _1 , and thus this 
source is classified at the near distance with quality grade B. 

24.78+0.08 - The recombination line velocity is 108.6 km 
s _1 , and Hi absorption is seen up to 108.0 km s _1 . Since the 



tangent point velocity is 112.8 km s 1 , we classify this source 
at the near distance with quality grade B. 

25.70+0.04 - There are two recombination line systems at 
this position, at 53.3 km s _1 and 102.0 km s _1 . The methanol 
maser emission ranges from 89 to 101 km s _1 , and is thus as- 
sociated with the 102.0 km s recombination line. The 21 cm 
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continuum image shows a single source with some sub-structure. 
Hi absorption at greater than 3<x level is seen up to 120.9 km 
s _1 , while the tangent point velocity is 1 12.4 km s _1 . This would 
place either recom bination line source at the far distance with 
quality grade A. iDownes et a l. (1980) list a single HllOo' ve- 
locity of 52.0 km s _I which suggests that the region seen at 21 
cm is not associated with the methanol maser. Thus, we do not 
give a distance to the methanol maser, and classify the 53.3 km 
s _1 component to be at the far distance. 

28.30-0.38 - L89 report two recombination lines at 44.8 km 
s _1 and 75.6 km s _I . The methanol maser emission, ranging from 
79 to 93 km s _1 is thus associated with the 75.6 km s _1 compo- 
nent. There are two distinct continuum sources (labeled A and B 
in Figure 1) seen in the 21 cm image, and both display H i spec- 
tra that place the source at the far distance. Hence, while we are 
unable to determine which continuum source gives rise to each 
recombination line, we classify both sources at the far distance. 

30. 70-0.07, 30. 76-0.05 and 30.82-0.05 - The three regions 
are part of W43, a large and complex region as seen in Figure 
1 . Absorption features are seen at velocities greater than the tan- 
gent point velocity of 104.5 km s 1 for 30.70-0.07 and 30.82- 
0.05, while for 30.76-0.05, absorption is seen up to 101.5 km 
s _1 . The systemic velocities for 30.70-0.07, 30.76-0.05 and 
30.82-0.05 are 101.0, 91.6 and 105.7 km s" 1 respectively. A 
straightforward application of the KDA resolution would put 
30.70-0.07 and 30.76-0.05 the far distance (with quality grades 
A and B respectively), and 30.82-0.05 at the tangent point. A 
stronger source in the same complex, 30.78-0.03 is also found 
to be at the far dista n ce (w ith a systemic velocity of 91.6 km 
s 1 ) by KolpaketaLl d2003h . Since the regions are part of the 
same complex, their different velocities can be explained as aris- 
ing from the random motions within the complex. Hence, we 
assign all sources at the far distance corresponding to 97.5 km 
s which is the mean velocity of the different components. 

32.97+0.04 - The recombination line velocity measured at 
GBT is at 93.2 km s~', whi le the tangent point velocity is 99.0 
km s" 1 . IWatson et al.l d2003l) do not classify the nearby (and pre- 
sumably the same) source 32.99+0.04 since its systemic velocity 
is within 10 km s _1 of the tangent point velocity. Hi absorption 
can be seen up to at least 104.1 km s _1 . Thus, we classify this 
source at the far distance with quality grade B. 

33.13-0.0 9 - This source is classi fied as being at t he tan - 
gent point by lArava et al.1 (120021) and iKuchar & Banial (11994 . 
The L89 recombination line velocity is at 93.8 km s _1 , and Hi 
absorption is seen up to 102.8 km s _1 . Since the tangent point ve- 
locity is 98.6 km s _1 , we classify this source at the far distance 
with quality grade B. 

33.40+0.01 - The GBT spectrum shows a prominent line 
at 72.3 km s _1 , and a weaker component at 105.4 km s _1 . The 
methanol maser, whose emission ranges from 96 to 107 km 
s _1 is associated with the latter component. The 21 cm contin- 
uum shows a compact source with a much weaker component to 
the southwest, the latter being too weak to extract a H i spectrum. 
The tangent point velocity at this Galactic longitude is 97.9 km 
s" 1 , which would place the methanol maser at the tangent point. 
IKuchar & Banial (U9 94) list the nearby (and potentially the same) 
source G33. 4-0.0 (with Galactic latitude and longitude of 33.42 
and 0.0 respectively) to be at a velocity of 76.5 km s and being 
at the far distance. It is thus possible that the brighter continuum 
source is responsible for the 72.3 km s recombination line, and 
the weaker recombination line arises from the weak continuum 
source or is undetected at 21 cm. Hence, we classify the 72.3 km 
s _1 component at the far distance, and locate the methanol maser 
at the tangent point. 



48.90-0.27 - This source, at a systemic velocity of 66.5 
km s _1 , is part of the W51 complex. Kinematically, this source 
would be at the tangent point since the tang ent point velocity is 
62. 2 km s" 1 (and has bee n classified so by |Kolpak et ail 2003 
and lKuchar & Banialfl994l) . However. iGenzel et all d 198 lb mea- 
sured a distance of 7 + 1 .5 kpc using statistical parallax by VLBI 
measurements of H2O masers in W51 MAIN. Since this is a 
measurement that is independent of the Galactic rotation curve 
or other distance indicators, we adopt this as the distance to this 
source. 

50.29+0 . 69 - The HI 10a recombination line measured by 
lArava et ail (|2002|) is at 25.9 km s _1 , while the methanol maser 
emission extends from 26 to 33 km s" 1 . This sour ce is classi- 
fied at the near distance by Watso rTet al.l (120031) and lArava et al.l 
(2002) using formaldehyde absorption. The H 1 absorption spec- 
trum however clearly places this source at the far distance with 
quality grade A. This is an example (albeit rare) of how the lack 
of formaldehyde clouds between the near and far distance points 
can give rise to incorrect classification of a source. 

81.87+0.78 - The methanol maser emission in this source 
ranges from to 13 km s _1 while the recombination line velocity 
is 1.7 km s _1 . The Hi spectrum shows absorption up to 15.2 
km s _1 , and the tangent point velocity is 3.0 km s _1 . Hence, we 
classify this source at the far distance with quality grade B. The 
use of the modified technique to calculate kinematic distances, 
described in this section 4, would however place this source at 
the tangent point, as is noted in Table 1 . 

6. Discussion 

16 out of 38 sources in our study with KDA resolutions have 
been targeted in previous studies. Sources for which our classifi- 
cation disagrees with that of previous work are noted in section 
5.1. Figure 2 shows the locations of the sources in a face-on dia- 
gram of the Galaxy. We als o compiled a list of H 11 r e gions whose 



KDA has been resol ved bylKuchar & Banial d 19941) . Arava et al 
(2002}, iFish et all d^OlLlKolpak et all (120031) IWatson et al 



(120031) . and lSewilo et al.l d2004 in addition to those in our study. 
All the distances were recalculated using the C85 rotation curve, 
and the modified technique described in section 4. Some KDA 
resolutions, especially of sources with velocities near the tan- 
gent point velocity, were reinterpreted by the criteria described 
in section 4 for consistency. Figure 3 shows all the sources in this 
compilation (we do not distinguish the sources by author since 
a number of sources are shared by multiple groups) superposed 
with two popular spiral ar m models of the Galaxy - the NE200 1 
model oflCordes & Laziol d2002l) andthemodel o flValled d2002l) ; 
IValled d 19951) . Figure 4 shows the same sources, but with dis- 
tances calculated using the modified procedure. The NE2001 
model defines spiral arms as loga rithmic spirals with pert urba- 
tions similar to those described in Tayl or & Cordesl £[993). The 
Vallee] d2002l) ; rValle e (1995) model defines the spiral arms as 
logarithmic spirals, the parameters of which are determined from 
a meta-an alysis o f nume r ous observatio nal datasets. It is to be 
noted that IValled d2002l) ; IValled (1 19951) use an R of 7.2 kpc, 
and consequently, we have scaled the model to Ro of 8.5 kpc for 
comparison with our data and the NE2001 model. 

The data points in Figs. 3 and 4 sample the Crux-Scutum, 
Carina-Sagittarius and Perseus arms. Qualitatively, the NE2001 
model is a better descrip tion of the G alactic spiral structure com- 
pared to the model of [ValleeJ (120021) . The discrepancy between 
the two models is greatest for the Carina-Sagittarius arm. A 
number of H 11 regions also appear to lie in inter-arm regions, 
although the error bars in the distances are almost always large 
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Fig. 2. Positions of the 6.7 GHz methanol masers and associated 
H ii regions for which we have resolved the KDA. The left panel 
shows the distances using the conventional method labeled d in 
Table 1, while the right panel shows the distances calculated us- 
ing the modified LSR velocity and rotation curve, labeled as d' 
in Table 1 . The open and filled triangles show the distances with 
quality grades A and B respectively. The location of the Galactic 
Center is marked as "GC". 



enough for these regions to be potentially associated with one or 
two arms. The lack of data points along the Outer arm could be 
a sensitivity effect, as only the brightest H n region complexes 
would have radio continuum high enough to be targeted in these 
stu dies. 

iFish et ail d2003l) proposed using the small vertical height of 
compact H n regions to resolve the distance ambiguity using a 
simple relation. Based on a Monte Carlo simulation, they sug- 
gest that H ii regions can be classified to be at the near or far dis- 
tance based on which distance was closer to a predicted value, 
dp = 1.84 \b\ . For our sample, this method has an accuracy 
of 70%. Surprisingly, the accuracy drops to 60% when restricted 
to sources that h ave Ad = df ar - d near > 4.7 kpc. In contrast, 
IFish et ail d2003l) found the method accurate to better than 90% 
for a sample with similar criteria in their simulation. It is to be 
cautioned that our sample has numerous biases, so that it can- 
not b e used to test the va lidity of the method, although we note 
that IWatson et al.l (|2003) had a similar success rate with their 
sample. Nevertheless, it indicates the danger of using this tech- 
nique for KDA resolution of specific sources as opposed to a 
statistical study of a large sample. As su c h, this specific line of 
argument used by ICorbel & Eikenberrvl d2004l) to justify their 
placing 10.32-0.15 at the far distance and 10.62-0.38 at the near 
distance is invalid. 

It is also curious that we find about 70% of the sources to be 
at the far distance. Th is is somewhat similar to the observation 
of lKolpak etai]j2003b . However, our methanol maser sample is 
not flux limited, and has a numerous selection effects and biases 
associated with it. Hence, we do not attempt to draw any conclu- 
sion from this observation. 



6.1. Application to Fainter Sources 

As mentioned earlier, the sources selected for inclusion in our 
study involved methanol masers in the first Galactic quadrant 
and at / > 10°, with direct or nearby associations of NVSS 



sources that were stronger than 30 mJy. This resulted in a se- 
lection of fewer than 30% of the masers in this region from the 
General Catalog, of which about 40% were too resolved in our 
data to solve the KDA. To extend this study to a larger sample 
of methanol masers, one would need to observe fainter sources. 
This raises the question of plausibility of applying the technique 
to faint continuum sources. 

One faces two challenges when dealing with faint sources. 
Firs, there is significant contamination from the stronger sources 
(which are numerous in the Galactic Plane) through sidelobes 
of the dirty beam. Second, there is low-level H i emission in the 
entire image. As explained in section 3, CLEANing the data is 
essential to be able to remove the contamination from stronger 
sources. However, a full field CLEAN does not work well for 
faint sources that are significantly contaminated by stronger 
sources in the field. On the other hand, a CLEAN with boxes 
at the locations of the continuum emission will tend to bring H i 
emission, which has a distribution that is completely different 
from that of the continuum, into the CLEAN boxes and into the 
target source. 

The problem of emission can be mitigated to a certain extent 
by going to a larger array configuration, which tends to filter 
out more of the extended Galactic emission. Doing this would, 
however, resolve the sources even further, requiring significantly 
greater telescope time to achieve sufficient signal to noise. This 
is especially so, since a 3cr detection of absorption at optical 
depths corresponding to exp(-r) of at least 0.3 is required to 
resolve the KDA with high confidence. A lower configuration 
(D or the proposed E configuration) would yield higher signal 
to noise ratio on the sources by keeping them unresolved, at the 
expense of introducing a larger extent of H i emission. 

Extending H i absorption against radio continuum emission 
for resolution of the KDA is thus a significant challenge in the 
case of weak continuum emission. A promising variant is to 
employ H i self-absorption to distinguish between near and far 
distances, a technique which has recently been use d to resolve 
the KDA towards MYSOs by iBusfield et alj d2006l) . Since this 
method relies on the absorption (or the lack thereof) of back- 
ground H i line emission at the systemic velocity of the source, 
it does not suffer from the problem of lack of 21 cm contin- 
uum. The only challenge is in the interpretation of the H i fea- 
tures and distinguishing self-absorption from the often complex 
emission line profile. Busfi eTdet al.l d2006l) estimate a success 
rate of ~ 70% for KDA resolution using this technique for their 
sample of MYSOs selected from the red MSX source (RMS) 
survey. Formaldehyde absorption at 6 cm may also be viable 
for sources having weak 21 cm emission. An advantage is that 



young H n regions, whose flux density has a 



dependence 



on frequency would be brighter at 6 cm, but a potential problem 
for the weaker sources is the increased risk of confusion from 
absorption against the CMB. 

7. Conclusion 

We have resolved the KDA towards 34 6.7 GHz methanol masers 
and 4 additional H n regions using H i absorption measured with 
the VLA. Over 70% of our sources have been observed to be at 
the far distance. T he use of Galactic latitude to resolve the KDA 
(IFish et alj|2003l) is not very successful for our sample. The es- 
timates of distances to the methanol maser sources obtained by 
resolving the KDA are crucial for the ongoing project to deter- 
mine the methanol maser luminosity function. 
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Fig. 1. Images of the 21 cm continuum, Hi optical depth spectra and GBT recombination line spectra (where applicable) of the 
sources in our study. The left panel shows the continuum image with contours at -2 and 2 times the rms noise listed in column (9) of 
Table 1, and at factors of 2 thereafter. For sources with multiple velocity components, the component from which the Hi spectrum 
is derived is marked with an 'X'. The right panel shows the line to continuum ratio, or exp(-r). The source velocity, v s , and the 
tangent point velocity, v t , are indicated by short and long dashed lines, respectively. For sources with GBT recombination line data, 
the top right panel shows the H i spectrum, while the bottom right panel shows the recombination line spectrum in emission. 
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X (kpc) X (kpc) 

Fig. 3. H ii regions wi th resolved KDA compiled from lKuchar & Banial (ll994l) . lAraya et alJ j2002h . lFish et all (l2003l) . lKoipak et al.l 
(2003), Watson etaD (120031) . ISewilo et alJ (2004) and our study. Distance s have been calculated using the C85 rotation curve. 
The left and right panels show the spiral arm loci of the NE2001 model of Cordes &Laziol dl002), and the Vallel d 19951) model 
respectively. For clarity, no error bars are shown. 
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X (kpc) X (kpc) 

Fig. 4. H ii regions with reso l ved KDA compiled from fKuchar & Banial (ll994h . lAraya et alJ j2002h . lFish et all (l2003h . lKorpak et al.l 
d2003lh IWatson et al.l d2003l) . ISewilo et alJ y004) and our study. Modified kinematic distances have been calculated using the pro - 
cedure described in S ection 4. The left and right panels show the spiral arm loci of the NE2001 model of Cordes & Laziol (|2002), 
and the lValleel {l995) model respectively. For clarity, no error bars are shown. 



